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ABSTRACT 

Context. The correlation between Ca and Ha chromospheric emission, known to be positive in the solar case, has been 
found to vary between -1 and 1 for other stars. 

Aims. Our objective is to understand the factors influencing this correlation in the solar case, and then to extrapolate 
our interpretation to other stars. 

Methods. We characterize the correlation between both types of emission in the solar case for different time scales. Then 
we determine the filling factors due to plages and filaments, and reconstruct the Ca and Ha emission to test different 
physical conditions in terms of plage and filament contrasts. 

Results. We have been able to precisely determine the correlation in the solar case as a function of the cycle phase. We 
interpret the results as reflecting the balance between the emission in plages and the absorption in filaments. We found 
that correlations close to zero or slightly negative can be obtained when considering the same spatio-temporal distri- 
bution of plages and filaments than on the sun but with greater contrast. However, with that assumption, correlations 
close to -1 cannot be obtained for example. Stars with a very low Ha contrast in plages and filaments well correlated 
with plages could produce a correlation close to -1. 

Conclusions. This study opens new ways to study stellar activity, and provides a new diagnosis that will ultimately help 
to understand the magnetic configuration of stars other than the sun. 
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1. Introduction 
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For the sun, the correlation between the emission in the Ca 
K and H lines and the e mission in Ha is wellkno wn, as both 
follow the solar cycle ([Livingston et al.l . |2007[ ). However, 
this may not be the c ase for other stars, as shown by 
ICincunegui et al.l (|2007[ ). who found that there was a large 
dispersion of such correlations. They have studied 109 stars, 
with a number of observations for each star ranging from a 
few to a maximum of 22 and covering about 7 years (total 
span of their observations). They found that for some stars, 
the correlation could be close to and even negative. They 
do not provide any interpretation of this result. 

There is therefore a need to understand the nature of 
this correlation. This type of observation can indeed be 
used as a diagnosis of stellar activity, but up to now no 
study has attempted to describe what kind of diagnosis 
could be derived. What can this correlation teach us about 
the nature of the stellar activity or the type of dynamo, for 
example ? 

To understand this, in this first paper (future work will 
be devoted to more accurate determination of this correla- 
tion for a large sample of stars), we study in detail this cor- 
relation in the case of the sun. We will first consider its de- 
pendence on cycle phase and the influence of the time span 
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used to compute it. The ICincunegui et al.l (|2007l ) observa- 
tions covered a range smaller than the solar cycle. When 
considering a short time span, we expect that a smaller 
correlation will be measured because of the small range of 
variation of the emission compared to the variability over 
short time scales . We u se the measurements performed by 
[Livingston et al.l (|2007f ) at Kitt Peak and covering 19 years 
to perform this study. This will allow us to determine if 
such a correlation varies over the solar cycle. 

Then we attempt to understand this correlation by the 
contribution of emission in plages in Ca and Ha, which 
should lead to a correlation very close to 1, and of absorp- 
tion in filaments, which should decrease the correlation be- 
cause many of them are not associated directly with active 
regions and because their contrast is much lower in Ca than 
in Ha. We use Meudon spectroheliograms obtained between 
1990 and 2002 in both Ca and Ha in order to estimate the 
relative contribution of plages and filaments in the solar 
case, and then to study their influence on the correlation 
between the two emissions. This provides clues to study the 
physics of stellar cycles, as it gives access to more informa- 
tion than the average activity level. 

The outline of this paper is as follows. In Sect. 2, we 
study in detail the correlation between the Ca and Ha 
emission in the solar case, using Kitt Peak observations. 
In Sect. 3 we study Meudon spectroheliograms to estimate 
the contribution of plages and fllamcnts to the observed cor- 
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relation. These results are used in Sect. 4 to in terpret the 
stellar observations of ICincunegui et al.1 (|2007D . We con- 
clude in Sect. 5. 



2. Analysis of Kitt Peak observations 

In this section, we consider the sun as a star in order to 
study the correlation between the Ca and Ha emissions. For 
this purpose, we use the long time series of observations per- 
forme d at Kitt Peak by W. Livingston (Living ston et al.l 
[2007lfl . The sun was observed as a star (full-disk integrated 
light) using the Fourier Transform Spectrometer at Kitt 
Peak, thus providing spectra with a very high spectral res- 
olution (about 500000). The intensity in the line core of 
several lines is also available on the web^. 



2.1. Ca and Ha emission variations 

We consider the intensity in the core for each of 
these lines, normalized by the continuum intensity (see 
[Livingston et al.l . l2007l for more detail). Measurements 
cover the range 1984-2003 and consist of 383 observations. 
They therefore cover one and a half cycles, and are evenly 
spread over the full time range. Fig. [1] shows the emission 
in both lines as a function of time. These variations are 
well correlated for example with the activity level defined 
by th e spot number, as already known (Livingston et all . 
|2007() : the correlation with the daily spot number is 0.86 
for Ca and 0.78 for Ha. 



2.2. Correlation between the Ca and Ha emissions 

We now study in more detail the relation between these 
emissions. We define the correlation factor between the two 
emissions Eca and Eho as : 
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We find a correlation of 0.806. The slope of the Ca emis- 
sion versus the Ha emission is 0.27. This means that there 
is a strong correlation between the two, but it is different 
from 1. Departure from 1 could be due to noise in the data 
and/or to physical causes, which will be discussed later. 

In Fig. [21 we study how the correlation varies with the 
phase of the solar cycle. We compute the correlation be- 
tween the two emissions for different time lags and different 
starting times (i.e. different cycle phases) in the series. For 
each time span, all dots correspond to a different starting 
time in the series, covering the whole range (minus the time 
span). For short time spans, very small correlations can 
be observed, including correlations close to zero or slightly 
smaller (^-0.3), i.e. correlations are biased toward values 
smaller than the one obtained on long time scales. The cor- 
relation for our whole time serie is reached only for time 
spans close to the cycle length or larger. 

Fig. [3] shows the variation of the cor relation over time , 
for a time span of 7 years (as used by ICincunegui et all . 
12001 . There is a strong relation between this correlation 
and the cycle phase. The correlation is indeed higher than 
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Fig. 1. Upper panel : Emission in the Ca line versus time 
(from W. Livinsgton). Lower panel : Emission in the Ha 
line versus time (from W. Livingston). 
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Fig. 2. Emission in the Ha line versus emission in the Ca 
line (from W. Livingston), from the 1984-2003 complete 
series. The correlation is 0.806 and the slope 0.27. 
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Fig. 3. Correlation between the emissions in the Ca and 
Ha lines, versus the time span. Each dot corresponds to a 
different cycle phase. The solid line represents the average 
and the dashed line the median. The horizontal dotted line 
corresponds to the correlation computed on the whole data 
set. 



average at the end of the ascending phase of the cycle 
(above 0.9) and much lower (about 0.45) at cycle minimum. 
Furthermore, the curve is not symmetric. For a given ac- 
tivity level, the correlation is higher during the ascending 
phase of the solar cycle compared to the descending phase. 
The behavior is similar for the slope. 
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3. Meudon spectroheliogram analysis 

3.1. Contributions to the correlation between the Ca and Ha 
emission 

We now explore possible causes for a correlation of less 
than 1 between the two emissions in the solar case. Active 
regions both show a strong emission in Ca and Ha due to 
plages and their position on the solar disk is similar, as 
shown in Fig. [5] However, the surface covered by the Ha 
plages is slightly smaller than the surface covered in Ca, 
as the latter are formed higher in the chromosphere. The 
chromospheric network is also more prominent in Ca than 
in Ha. We therefore expect a small dispersion between the 
two filling factors, that should lead to a correlation slightly 
lower than 1. 

An important difference between the two is the presence 
of dark filaments at the surface, mostly visible in Ha. They 
are also slightly absorbing in Ca, but with a much lower 
contrast. When considering the integrated emission, their 
contribution will therefore be significant for Ha but not for 
Ca : they will decrease the correlation between the Ca and 
Ha emissions towards zero, because they do not coincide 
exactly with plages. On the other hand, if they were well 
correlated with plages, a strong filament contrast would 
lead to an anti-correlation. Filaments are cool and dense 
gas structures in the corona, maintained against gravity by 
the presence of aa magnetic field. They are formed at the 
boundary between regions of opposite-polarity line-of-sight 
mag netic fields, ei ther in the quiet sun or in active regions 
(see iMartinl . Il998l . for a review). 
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Fig. 4. Upper panel : original emission in the Ca line versus 
time (stars) and average over a time span of 7 years (solid 
line). The time corresponds to the middle of the 7 year 
span. Middle panel : Correlation between the Ca and Ha 
emission versus time for a time span of 7 years. The hori- 
zontal dashed line corresponds to the full data set. Lower 
panel : Slope of the Ha emission versus the Ca emission for 
the same conditions. 



There are two different categories of filaments. Some 
are present in active regions, and these should naturally 
correlate well with the activity cycle both on short (days) 
and long time scales (cycle). However, they are usually 
small and very active, and therefore do not live very long. 
Therefore their contribution to the correlation should be 
small, as these properties are not compensated by very large 
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Fig. 5. Upper panel: A Ca spectroheliogram (left) and the corresponding segmented image (right). The outer ring shows 
the solar limb, and the inner ring the limit at 0.9 xR© used to select structures (plages in that case). Lower panel: Same 
for Ha, with the bright structures corresponding to plages and the grey structures to filaments. 



numbers. Furthermore, those that are stable in active re- 
gions represent a small number, and they seem to show n o 
correlation at all with the solar cycle (jMackav et al.l . [2008( 1. 



Secondly, there are also many large and stable qui- 
escent filaments outside active regions. These filaments 
should therefore contribute significantly to the corre- 
lation. They s eem to be corr e lated with the solar 
cycle dMouradian fc Soru-Escautl 11994 iLi et all . l2007l 
iMackavet al.l . 12008( 1 . but again at a large scale (the cycle) 
and not necessarily at the scale of days or weeks. Even if 
their presence is related to active regions and in some way 
ultimately to large regions that are predominantly unipolar, 
the correlation may be low on small time scales, as there is 
a time delay between the active regions and the formation 
of these unipolar network regions. 



3.2. Data analysis 

To estimate the contribution of plages and filaments, we 
measure their filling factor over the solar cycle using 
Meudon spectroheliograms (provided by BASS200CQ). As 
mentioned above, the filament contrast in Ca is small, and 
we neglect it in the following. We therefore measure the 
plage filling factor both in Ca (taken at the center of the 
Ca II K line and showing the upper layer of the chromo- 
sphere) and in Ha images, and the filament filling factor 
in Ha images only. We have used numerized data obtained 
between 1990 and 2002, i.e. covering an 11-year solar cycle. 
Fig. [S] shows the number of image pairs (Ca and Ha) for 
each 30 day period. There are 1690 pairs in our data set. 
There are naturally more images during the summer, due 
to the better weather. There are no large gaps in the data 
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Fig. 6. Number of image pairs in the Meudon data set 



set, except for a few months in 1996, and the whole period 
is therefore well sampled. 

It is beyond the scope of this paper to perform a very 
sophisticate d regio n extraction s uch a s those realized by 
[Fuller et alJ (l2005t). Ilpson et alj (I2005D . iBernasconi et al " 



(120051 ) ■ lAboudarham et all (|2008D or IScholl fc Habbal 



20081 ). as their purpose was to identify complete filaments 



as a whole, with Virtual Observatory applications in mind 
(building of synoptic maps for example, in which the differ- 
ent pieces constituting a filament must be associated and 
tracked in time). Here we are mainly interested in the total 
area covered by these structures at a given time, therefore 
there is no need to associate them. This considerably sim- 
plifies the segmentation step. For the same reason, we do 
not attempt to use all data, and have therefore selected the 
best images (reasonable seeing, no clouds, no large dust 
line). Therefore our algorithm, although presenting strong 
similarities with these previous works, is more simple. It 
is constituted of the following steps: center-to-limb darken- 
ing correction; analysis of the intensity histograms in order 
to determine a threshold (ensuring robustness of the final 
filling factor over time); application of this threshold to 
the data; application of a thinning procedure; analysis of 
the segmented data to extract the size of the structures. 
Structures smaller than 10 pixels are eliminated (see fur- 
ther for more discussion on the size threshold). 

A size selection is performed to eliminate the smallest 
structures, as we focus on plages and active network to re- 
duce the noise. This size threshold is applied only to the Ca 
structures (and fixed at 120 Mm^). We have tested thresh- 
olds between 30 Mm^ (corresponding to about 10 pixels 
and therefore very noisy) and 270 Mm^ (corresponding to 
small active regions). The chosen threshold corresponds to 
the smallest structures that did not significantly influence 
the results and allows us to take into account a signifi- 
cant proportion of the magnetic network. As the very small 
structures have a lower contrast (jWorden et al.1 . Il998f ) , it is 
not a problem to eliminate them at this stage. Then we only 
keep structures that are overlapping in both wavelengths. 
This step makes the procedure robust by eliminating small 
features that would appear in one of them only due to the 
threshold choice, if their size were to be close to the thresh- 



old. For the same reason, we have selected filaments with 
the same threshold (sizes larger than 120 Mm^), but the 
infiuence is not significant either because the smallest fea- 
tures do have a very low contrast as well and therefore do 
not influence the integrated emission. As shown in Fig. [SI 
large filaments are well defined, however we do miss the 
thinnest filaments, which are present mostly in active re- 
gions. Because these are much thinner and short-lived, they 
are the filaments which contribute the less, so it is not a 
critical point (see Sect. 3.1 for a discussion). 

Finally, a filling factor is computed using the selected 
structures (surface divided by the total surface of the sun 
on the image) . We are mostly interested in the variability of 
these filling factors and the correlation between them, more 
than in the absolute values of the filling factors themselves, 
so the exact thresholds used to extract the structures is 
not as sensitive as it would be to get the exact size of these 
structures for example. 

3.3. Filling factors over the cycle for active regions and 
filaments 

Fig. [7] shows the resulting filling factor versus time for the 
3 types of structures: plages identified on Ca II K images 
(ffca), plages identified on Ha images (ffna), and filaments 
identified on Ha images as well (fffii). As expected, we ob- 
serve a strong correlation of 0.92 between the plage filling 
factor in Ca images and in Ha images. The filament filling 
factor is also higher at cycle maximum, but the correlation 
with the plage filling factor is lower, with values of 0.52 for 
Ca and 0.53 for Ha. 

Fig. [5] also shows Ahq and fffii versus ffca- The correla- 
tion is quite clear in both cases. For Ha plages, the relation 
is close to linear with a change in slope for ffca around 
0.04: at higher activity levels, Ahq increases faster than 
ffca- This may be due to the largest structures at that time, 
with therefore a ratio closer to one between the size of the 
regions observed in Ca and Ha. As for filaments, a linear 
increase is observed for ffca up to 0.06, and then a satu- 
ration (corresponding to cycle maximum): above a certain 
activity level, the filament filling factor remains constant, 
which is a new result. This and the large dispersion explain 
the correlation close to 0.5, i.e. far from 1. We also note that 
at cycle minimum, when ffca and Ahq both go very close to 
zero, fffii is significantly above zero, showing the presence 
of a significant number of filaments, as already observed 
(iMakarov fc Mikhailutsal . Il992l . iMouradian fc Soru-Escaul , 
I1994D . 

The different time scales will be important when inter- 
preting stellar data. On a long time scale (cycle period for 
example), the filament filling factor is roughly correlated 
with the activity level (determined from the Ca filling fac- 
tor fo r example), as already pointed out by iMackav et al.l 
(|20080 for a few selected periods over the cycle. In addi- 
tion to the saturation at high activity levels found above, 
the correlation is not close to 1 due to a large dispersion, 
and this means that at small time scales (rotation period) 
the correlation could be much closer to zero. Finally, Fig. [5] 
shows the filling factors for Ca and filaments for an aver- 
aging over temporal windows of about 4 rotation periods. 
On this plot, fffii has been normalized to ffca (same me- 
dian) to make the comparison easier. The correlation at 
this large temporal scale is now much closer to 1 (0.87) due 
to the temporal smoothing. The value of 0.92 (as for ffca 
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Fig. 7. Upper panel : Filling factor of Ca II K plages. 
Middle panel : Filling factor of Ha plages. Lower panel : 
Filling factor of filaments. 



and ffjjQ obtained above) is found for a temporal smoothing 
of about 200 days. There are, however, a few noticeable dif- 
ferences. At cycle minimum and during the beginning of the 
ascending phase, the surface covered by filaments increases 
faster than the plage filling factor. On the other hand, at 
cycle maximum, the peaks are not as sharp, especially for 
the 1990-1991 maximum, which is related to the saturation 
above. The cycle maximum in 2000-2002 is known to ex- 
hibit two peaks, as shown by the ffca : this is not observed 
for filaments, with in fact a peak in between, shifted by 
about 200-300 days from the first peak. Fig. [9] also shows 
that the dip in the ratio between the plage and filament 
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Fig. 8. Upper panel : Filling factor of Ha plages versus 
filling factor of Ca II K plages. Lower panel : Filling factor 
of filaments versus filling factor of Ca II K plages. 

filling factors observed around 1996-1997 corresponds to 
the low correlation factor observed in the Kitt Peak emis- 
sion at that time, which reinforces our interpretation of the 
correlation evolution by the presence of filaments. 

3.4. Reconstruction of the emission 

In using contrasts for plages in Ca and Ha, and for filaments 
in Ha (Cca7 Cho and Cfii respectively), it is possible to 
reconstruct the emissions in Ca and Ha as defined by : 

Eca = Cca X ffca 

Ehq = Cna X Shq — Cfil X fffii. 

In the following, we have assumed constant contrast over 
time. We have also considered the same contrast for struc- 
tures of all size. We know that for bright structures, there 
is a variabilit y depending on their size from plages down to 
the network ( Worden et al.l . [T998[ ) . however in this work we 
focus mostly on plages (see Fig. [S]), so we do not expect a 
large effect. 

Unfortunately, the available Meudon data alone do not 
allow us to compute the plage and filament contrast, as 
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Fig. 9. Upper panel : Ca filling factor ffca (solid line) and 
filament filling factor fffii normalized to the same median 
(dashed line), both averaged of 4 rotation periods, versus 
time. Lower panel : Correlation between the Ca and Ha Kitt 
Peak emission versus time (solid line, from Fig. 4, middle 
panel) superimposed on the ratio (dashed line) between the 
two upper curve (ffca/fffii with a normalization of fFfii as in 
the upper panel). 



Table 1. Correlations and parameters from a fit in the Kitt 
Peak data or from the literature. 



Type 


Variable 


Kitt Peak 


Fit 


Literature 


Obs. 


Corr(ffca,fiHa) 
Uorr^n:ca,n:fii ) 
Corr(ffHa,fffii) 




0.923 

U.0Z4 

0.520 


0.923 

U.0Z4 

0.520 


Par am. 


Cca 




0.20 


0.70 








0.13 


0.35 




Cfli 




0.30 


0.20 


Constr. 


Corr(Eca,EHa) 


0.805 


0.805 


0.913 




Slope(Eca,EHa) 


0.267 


0.268 


0.260 




Ampl. Eca 


0.028 


0.028 


0.096 




Ampl. Ehq 


0.012 


0.013 


0.027 








0.0009 


0.109 



Note 1 . The first 3 lines provide the correlation between the ob- 
served filling factors. The following set of lines gives the contrasts 
which have been used (either from a fit on from the literature). 
The last set of lines shows the constrains that can be used to 
either constrain the fit or, in the last column, to test the con- 
trasts found in the literature, compared to the values expected 
from the Kitt Peak observations. 
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Fig. 10. Reconstructed emission in Ca II K versus the re- 
constructed emission in Ha. 



there is no intensity calibration. It is, however, possible to 
determine them using both Meudon and Kitt Peak data. 
We have estimated the best contrasts Cca, Chq and Cgi 
that would lead to reconstructed emissions varying over a 
range similar to that observed at Kitt Peak. We also con- 
strain these contrasts by adjusting them in order to obtain a 
correlation and a slope between the reconstructed emissions 
as close as possible to the Kitt Peak values. We therefore 
minimize a which takes into account these three fac- 
tors, equally weighted. The smallest then provides the 
best estimates for the three contrasts. Using these contrasts 
and the observed filling factors, we derive the reconstructed 
emissions Eca and Eho- We find contrasts that reproduce 
the Kitt Peak observations (see Table [T]) . These contrasts 
correspond to a small band pass, as we are using Kitt Peak 
intensities in the line core. Fig. [10] shows Eh Q, versus E^a in 
that case and can be compared to the Kitt Peak observa- 



tions in Fig. [21 We note that despite the same amplitude, 
correlation and slope, the behavior is slightly different, as 
Kitt Peak observations exhibit the largest dispersion for 
small activity levels, while here the largest dispersion is for 
large activity levels. 

Another possibility is to use the contrasts published in 
the literature. Unfortunately, we did not find any large and 
systematic studies that would provide all necessary con- 
trasts. One reason is that they depend significantly on the 
band pass and that there are few studies on a large sam- 
ple. For C a structures, w e have used the value of 0.7 de- 
duced from lWorden et al] (|1998D . who obtained a variation 
with the type of str ucture (from the network to plages). 
lOrtiz &: Rastj (|2005f ) also observed a variation of the Ca 
contrast as a function of the line-of-sight magnetic field. 
We considered the value of .35 for Chq, assuming a facto r 
of two between Ca and Ha (jKononovich fc Nikulinl . [i975[ ). 
This ratio of 2 is close to what we obtained above with 
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Fig. 11. Comparison of the distribution of correlations ob- 
tained by Cincunegui et al. (2007) for a sample of 109 stars 
(solid line) and for the sun viewed from the equator nor- 
malized to the same number of stars (this work) for various 
phase of the solar cycle and a similar time span, i.e. 7 years 
(red line). The green line is derived from the same condi- 
tions as for the red line, but takes into account a uniform 
distribution of inclination angles between and 90° for a 
latitudinal extent of [-30°, 30°] for plages and [-60°, 60°] for 
filaments (see Sect. 4.1 for details). 

Table 2. Filament contrasts for which the correlation is 
equal to zero. 



Contrast 


Fit 


Literature 


Cca 


0.20 


0.70 


Chq 


0.13 


0.35 


Cfli 


0.30 


0.20 


Cfli(corr=0) 


1.16 


3.08 


Cm(corr=0)/Cfli 


3.9 


15.4 



Note 2. Given the fitted contrasts for the three types of struc- 
tures or those found in the literature, a filament contrast for 
which the correlation is equal to zero is given, in absolute value 
and relative to the solar value (see Sect. 4 for details). 



the fit. The filament contrasts show a very large disper- 
sion in the literat ure and we hav e ta ken an average valu e of 
0.2 deduced from lMaltbvl (| 19761) and lChae etall (|2006l ). As 
shown in Table [Tl the is much worse in that case, and 
it is not possible to reproduce the observed amplitude of 
variation and the correlation, as the plage contrasts are too 
high and the filament contrast too low. This may be due 
to the difference between the band pass at which the y were 
measured (0.5 A for example in the case of Wordcn et al.l . 
Il998() and the Kitt Peak observations we are using. 



4. Toward diagnosis of stellar activity 

4.1. Distribution of stellar and solar correlations 

Fig. [TT] sho ws the distri bution of correlations obtained by 
ICincuncgui et al.1 ()2007f ). They observe a significant num- 
ber of stars with correlations close to zero or negative. 
Superimposed on it, we show the distribution in the case of 



the sun, for a time span of 7 years (red line), and covering 
different phases of the solar cycle, deduced from the dots of 
Fig. [3] for that time span. Several conclusions can be drawn 
from this comparison. 

First, the two distributions are very different, as the tail 
toward small correlations is of small amplitude in the case 
of the su n, c ompare d to the nuin ber of cases around 0.8. 
However, [cinc uneg ui et al.1 (|2007l ) provide no uncertainties 
on the correlations, and some of these are large (given the 
available plots). However, large uncertainties probably do 
not suffice to explain the difference in distribution. This 
means that the activity of the sun is probably not typical 
of that of other stars, possibly due to filaments. A pos- 
sible source of discre pancy could be th e hete rogeneity in 
the sample studied bv lCincunegui et all (|2007f) in terms of 
stellar type. 

It is also necessary to take into account the large range 
of inclination angles present in stellar observations. In the 
solar case, the latitudinal extent of plages and filaments is 
different, so one can expect a different contribution from 
these two types of structures depending on the inclination 
angle. If considering that plages extend up to 30° and fila- 
ments up to 60°, for example, and assuming that the filling 
factor at each time step is attributed to that latitudinal 
range, a correlation of 0.8 (sun viewed from the equator) 
would fall to 0.75 for an inclination angle of 45° and to 0.49 
for a pole-on orientation. For a uniform distribution of in- 
clination angles between and 90°, the distribution shown 
in Fig. [11] (red line) should then be smeared toward lower 
correlations and is shown in green. 

The second conclusion is that despite these differences, 
the solar distribution exhibits small correlations in a few 
cases. This means that if stellar observations do not cover 
a full activity cycle (or more generally, the full range of ac- 
tivity existing for that star), we are likely to significantly 
underestimate the correlation with respect to the one that 
would be obtained when considering a full cycle. In that 
case, the only way to consider the data is a statistical anal- 
ysis of a large number of stars, i.e. via the distribution of 
correlations. Indeed, for a given star, the correlation may 
be low (for example close to zero) even if in reality the star 
has a correlation computed over the whole cycle similar to 
the solar one. If the cycle length is known, as well as the 
phase, it may be possible, in principle, to apply a correc- 
tion deduced from the plot of Fig. |4| (or a similar one if 
the time span of the observations is different). In that case, 
a distribution of corrected correlations could be produced, 
that would give an indication of the true correlations over 
the whole activity range. However, such corrections would 
assume a behavior similar to the sun (in particular, relative 
contrast). 

4.2. Influence of the filament contrast on the correlation 

We now consider cases for which the correlation is com- 
puted over the whole cycle. It would also be interesting to 
know what the sensitivity is of the correlation for filament 
contrasts differing from the solar contrasts found above. 
Fig. [T2I shows the correlation between the reconstructed 
emission in Ca and Ha for various filament contrasts, while 
keeping the same contrasts for plages in Ca and Ha. For 
low contrast, the correlation is naturally close to 1. For in- 
creasing filament contrast, the correlation decreases, crosses 
zero for a certain contrast, and then saturates for values 
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Fig. 12. Upper panel: Correlation between the recon- 
structed emission in Ca II K and Ha for various filament 
contrasts (relative to the fitted solar value). Lower panel: 
Same versus the actual filament contrast, without the nor- 
malization to the solar value. 



around -0.4. Table [5] shows the contrast necessary to reach 
a correlation of zero as well as the ratio to the solar con- 
trast. A contrast 3.9 times the solar value of 0.30 (i.e. 1.16) 
and the same distribution in time and position of the fil- 
aments would lead to a correlation close to 0. The factor 
is, of course, the same when applied to the area covered 
by the filaments instead of their contrast. When using the 
contrast obtained from the literature, the filament contrast 
necessary to reach a zero correlation is much higher. A con- 
trast of 1.16 is larger than 1, which means that a correla- 
tion of zero for the same spatio-temporal organization of 
the filaments can be reached only if the filling factor of fil- 
aments is increased. We see that, for example, for a factor 
2 increase in contrast and in size, which is not very ex- 
treme (the observed contrast of individual solar filaments 
varies by a factor of 3-4) , at least such a zero correlation is 
reached. 

However, we cannot obtain a strong negative correlation 
below ~-0.4, and such values correspond to extreme condi- 
tions compared to the sun. Therefore, the spatio-temporal 



organization of filaments must be significantly different if 
the presence of filaments explains the Ca - Ha relation ob- 
served in other stars. For example, filaments that would 
be more strongly in phase with plages and with a con- 
trast dominating that of Ha plages could lead to an anti- 
correlation. 

4.3. Influence of the filament and plage contrasts on the 
correlation 

We now look at the influence of all contrasts on the corre- 
lation and the slope simultaneously. Fig. IT51 shows the cor- 
relatio n versus the slopes as observed bv lCincunegui et all 
()2007D superimposed with the expected correlations and 
slopes for solar filling factors multiplied by various values 
of the Ha plage contrast (between and 1) and of the fil- 
ament (between and 2, i.e. allowing a larger filling factor 
for the filaments). Most stars could then be explained by 
this interpretation, i.e. with the assumption that plages and 
filaments have the same spatio-temporal distribution than 
on the sun but different contrasts. The most notable excep- 
tions are a few stars with correlations very close to -1. There 
are also a few stars with correlation between 0.92 (the solar 
case with no filament) and 1 but the uncertainties on the 
correlation in these stellar observations are larger than this 
difference. 

This graph showing the correlation versus the slope al- 
lows us to explore different physical conditions that may be 
found on stars of different types (e.g. different temperature 
and density profiles in the chromosphere). The value of the 
correlation depends significantly on several factors. First, it 
depends on the ratio between the plage contrast and the fil- 
ament contrast, as shown in Sect. 4.2. It also varies with the 
ratio between the plage contrast in Ca and Ha. These lines 
are not formed at the same level in the chromosphere, and 
we could expect stars of different types to exhibit a ratio 
different from the solar value. For stars with a Ha contrast 
much smaller than the Ca one, the correlation falls rapidly 
to zero or to negative values. 

5. Conclusion 

The first result of this work is a precise characterization of 
the correlation between the Ca and Ha emission in the solar 
case : the value of the correlation, dependence on the time 
scale and cycle phase, and finally interpretation using fila- 
ments, which decrease the correlation due to plages. Also, 
we have found that the filament filling factor is correlated 
to the activity level (0.5) but with a large dispersion (due 
to the low correlation at short time scales) and a saturation 
effect at high activity levels, which is new. 

In addition, we can use these results to estimate what 
would be expected on stars with a similar spatio-temporal 
distribution of filaments and plages but different sizes or 
contrasts. We can expla in most of the results obtained by 
ICincunegui et all ()2007D , but not all of them. This may re- 
flect physical conditions different from solar, either in chro- 
mospheric properties or in the activity c ycle. However, the 
uncert ainties in the results obtained by ICincunegui et all 
(|2007D seem to be quite large and they could be responsible 
for the observed discrepancy with the solar case. There is 
therefore a need to study in more detail the precise relation- 
ship between the Ca and Ha emissions in a large star sam- 
ple. This is probably at least partially due to the poor tern- 
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Fig. 13. Correlation between the Ca and Ha emission versus the slope, for the Cincunegui et al. (2007) results (stars) 
and the sun (circle). A model corresponding to the solar Cfii and a variable Chq is superimposed (thick line), as well 
as 5 models (thin lines), each of them with a different Chi (from a low filament contrast for the upper curve to a high 
filament contrast for the lower curve). For each curve, Chq varies from O.Of to 1. Cca is constant (fitted solar value of 
0.20). 



poral sampling. This is also necessary because their sample 
is quite heterogeneous in stellar types. This is a work in 
progress. 

If we want to extrapolate to stellar cases with differ- 
ent magnetic configurations and consequently a different 
spatio-temporal distribution of filaments, we need to obtain 
a better estimate of the filament behavior over the solar cy- 
cle. We therefore plan to study this u sing the very large dat a 
set of MDI/SOHO magnetograms (jScherrer et al.l . Il995^ . 
This will allow us to establish a link to the magnetic con- 
figuration. 
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